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Abstract 

Gel electrophoresis allows to separate knotted DNA (nicked circu- 
lar) of equal length according to the knot type. At low electric fields, 
complex knots being more compact, drift faster than simpler knots. 
Recent experiments have shown that the drift velocity dependence on 
the knot type is inverted when changing from low to high electric fields. 
We present a computer simulation on a lattice of a closed, knotted, 
charged DNA chain drifting in an external electric field in a topolog- 
ically restricted medium. Using a simple Monte Carlo algorithm, the 
dependence of the electrophoretic migration of the DNA molecules on 
the type of knot and on the electric field intensity was investigated. 
The results are in qualitative agreement with electrophoretic experi- 
ments done under conditions of low and high electric fields: especially 
the inversion of the behavior from low to high electric field could be 
reproduced. The knot topology imposes on the problem the constrain 
of self-avoidance, which is the final cause of the observed behavior in 
strong electric field. 
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Gel electrophoresis of linear and circular DNA and its dynamics has 
been since long time a topic on which numerical simulations and analytical 

models have been applied |H[2||2||2||5I|£|[ZII2|- 

Most experimental and theoretical studies of the electrophoresis process 
deal with linear or circular DNA [UJ EB EH E] • But DNA comes also in 
knotted form. Various classes of enzymes (topoisomerases and site-specific 
recombination enzymes) produce different types of knots or catenanes by 
acting on circular DNA molecules |13l I14j . The analysis of these knots gives 
some information about the mechanisms by which these enzymes are in- 
volved in the proper functioning of chromosomes (see for example |15| ) and 
about DNA packing Being able to study which knots are produced by 
a given enzyme in prescribed conditions implies being able to perform some 
sort of "knot spectroscopy", which can be done for example by electron mi- 
croscopy, where knots are observed one by one. Yet, if large numbers of knots 
need to be classified, then some high throughput technique is needed. Such 
a technique is gel electrophoresis. Indeed, experimental work has shown 
a linear relationship between the distance of electrophoretic migration on 
agarose gel of different types of DNA knots (all with the same number of 
base pairs) and the average crossing number of the ideal geometrical repre- 
sentations of the corresponding knots (closely related to the complexity of 
the knot)[T8]. As a consequence, the type of a knot can be simply identified 
by measuring its position on the gel, without the need of electron microscopy 
experiments as required until recently. 

At low electric field the usual observation is that the more complex the 
knot is, the higher is its mobility. A simple intuitive explanation for this 
behavior is that the compactness of a knot increases with its complexity (for 
a constant string length) and the friction coefficient 7 = Gitt/qRh (with Rh 
the hydrodynamic radius of the knot and rj the viscosity of the solvent) 
is smaller, resulting in higher mobilities. A more refined calculation of the 
friction coefficient 7 relies on Kirkwood-Riseman formula |19j : 



where the chain is modelled by iV beads of radius Ui and friction coefficent 

= 67Tr] ai, and Rij is the distance between beads i and j. The term 
due to hydrodynamic interactions between beads, in the second factor of 
equation Q explains the observed behavior: more compact molecules have 
smaller distances Rij, and thus a smaller friction coefficient. The calcula- 
tion of an average friction coefficient 7 on an equilibrium set of thermally 
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agitated DNA molecules forming different types of knots has confirmed the 
experimental results |2"Uj . 

Recently, it was observed by two-dimensional agarose gel electrophoresis 
that when the strength of the electric field is increased, the electrophoretic 
mobility of DNA knots changes behavior (fig. U) |171 I21| . The experiment 
was performed in two steps: a low strength electric field of 0.6 Vcm~ 1 was 
first applied along one direction in the gel. At this step, more complex knots 
show a higher mobility, in agreement with Kirkwood-Riseman formula. The 
same procedure is repeated in a second step but with a stronger electric 
field (f>Vcm~ l ) applied perpendicularly to the first one. In this case, the 
opposite behavior is observed: more complex knots cover smaller distances 
than simple ones. 

The presence of two regimes of weak and strong electric field can be 
captured with a simple argument. A knot of size £ drifts over a distance equal 
to its size along the direction of the applied electric field £ in a time t% — 
£/v = j£/q£, where v = q£/ r y is the drift velocity in the stationary regime 
and q the total electrical charge carried by the DNA molecule. During the 
same time, the drifting knot diffuses laterally over a distance d ~ yj2Dt^ = 



Boltzmann's constant. If the transverse diffusion explores distances d much 
larger than the typical size of the knot, £, then, on average, the knot will be 
able to avoid a collision with the gel and the knots will drift as if they were in 
a pure solution, with just a slight renormalization of the friction coefficient. 
If instead d <C £, then whenever a knot is on a collision course with a gel 
strand, it can not avoid it. As a consequence, after impinging over the 
obstacle, the knot needs to crawl around it in order to free itself. Crawling 
around an obstacle is much more difficult for more complex knots than for 
simple ones, due to the self- avoidance constraint. Following this argument, 
the two electric field regimes are separated by a critical field £ c that can be 
obtained by setting d ~ £, giving £ c ~ 2^^. In order to estimate £ c we 
use the typical values for a ll'OOO base pairs DNA knot: the size of a closed 
DNA ring is about £ = 300 nm |22|: the total charge q depends on the gel 
conditions, since every nucleotide carries a P~ group, hence one electron 
negative charge, that can be strongly screened by charges in the solvent. As 
a consequence we use q ~ 10~ 15 — 10~ 16 C We then obtain a critical field 
£ c = 0.1 — 1 Vcm" 1 in reasonable agreement with experiments |171 l2"Tj . The 
expression for the critical electric field holds also in the case when the gel 
is concentrated as it is the case in many experiments. Under the condition 
of high gel density, the DNA is filling the pores and the expression for the 




where T is the absolute temperature and ks is the 
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collision condition is that the DNA lying between two gel strands will collide 
with one gel strand before the DNA can drift transversally to the electric 
field. Instead of £, one has to insert the gel pore size I. The lateral diffusion 
constant D* has to be rescaled in order to include the effect of the gel. The 
condition for the critical electric field reads again: £ c ~ 2^|p, 

One has also to note that the two conditions for the critical electric field 
£ c actually mean that the energy gained by the DNA when moving one 
diameter £ or one pore size £ along the electric field is equal to twice the 
thermal energy: q£ c £ (or £) ~ 2/c#T. Although the above model gives a 
first hint of the origin of the observed behavior, here we want to address the 
issue more thoroughly using lattice Monte Carlo simulations. 

DNA knots are modelled by closed self-avoiding walks (SAWs) composed 
of N segments of length a on a three-dimensional cubic lattice (the lattice 
constant a is comparable to the persistence length of the DNA molecules). 
The gel is a two dimensional grid forming a sublattice with a mesh size b 
(= gel parameter) and perpendicular to the applied electric field (so that 
no knots can ever get impaled). The gel lattice is shifted by the quantity 
(f ' f ' f ) compared to the knot lattice, so that no points of the knot lie on 
the gel. Knots are not allowed to cross the gel network. The coordinates of 
the N monomers in the configuration at time t are written as: 

r(t) = (f 1 (t),r 2 (t),...,f N (t)) (2) 

with constraints ||r}(i) — r} + i(i)|| = a. 

The dynamics is followed using the BFACF algorithm [23]. Two types 
of moves are allowed: (a) the creation/destruction of a handle and (b) the 
flip of a corner into the mirror position (see Fig|5J). The first move clearly 
does not preserve the knot length, which can vary by ±2 at every step, but 
introduces the knot elasticity. The BFACF algorithm preserves knot classes, 
within which it is ergodic [21]. Self-avoidance is imposed by disallowing 
monomers to visit any site which is already occupied by other monomers. 
Furthermore, knots are not allowed to cross gel rods, so that corner flips and 
handle creation/destruction are forbidden when a rod has to be crossed. 

Under an external uniform electric field £, the electrostatic energy at 
time t is given by: 

N 

s ff(*) = -]vfe ( 3 ) 

N(t) is the length of the knot at time t, and it is associated with an elastic 
energy 

E el (t)= l -K[N{t)-N Q \ 2 (4) 
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where K is the spring constant. In the simulation a value K/ksT = 0.1 was 
used. The knot energy is then E(t) = E q (t) + E e i(t). 

At each timestep, we choose a point at random on the chain and pro- 
pose alternatively one of the two moves. If it satisfies the self-avoiding 
and gel- avoiding constraints, it is accepted with a probability given by the 
Metropolis algorithm: if the energy of the new trial configuration, Et r i a i, 
is lower than that of the previous configuration, E Q id = E(t), the move is 
accepted and f(t + 1) = f trial', otherwise, the probability of acceptance of 
the trial configuration is equal to exp{—[E tr iai — E(t)]} /ksT '. If the move is 
rejected, then f(t + 1) = f(t). 

After a knotted configuration is randomly generated, the knot type is ob- 
tained by calculating its Alexander polynomial 25]. Then, we let the system 
freely relax to thermodynamic equilibrium in the absence of an external field 
(£ = 0) until correlations from the initial configuration have disappeared. 
Then the electric field is switched on, and we let the knots migrate on the 
lattice. The quantities we compute are the position of the center-of-mass 
and the average crossing number (ACN) of the knot along a trajectory. 

Time is measured in Monte Carlo iterations, length in lattice spacing. 
The initial length iVo of our polymers was set to 150, and the mean length 
of the knot depends generally on the electric field and on the gel parameter. 
However, the mean length is 145 (146) for C=0.1 (C=0.4) and b=20 and 
we checked that during the simulations it fluctuates around that value. The 
average length is slightly shorter that Nq, since the probability of shortening 
the polymer is a slightly larger than the probability of lengthening it due to 
the self-avoiding condition. The gel parameter was set to b = 5, 10, 20 (in 
units of a), corresponding to a relatively sparse gel with big pores. For each 
initial knot, 20 • 10 6 iterations were performed. The center-of-mass position 
has been measured every 1000 Monte Carlo steps, and it was then averaged 
over the trajectories obtained by the migration of 100-200 different initial 
knots (to obtain an accuracy of about 10%). 

One problem with the Monte Carlo algorithm is that more complex knots 
have a smaller drift velocity than less complex ones even in absence of the 
gel, when time is measured in Monte Carlo steps. This is due to the fact that 
already in the absence of the gel, the moves are hindered by the complexity 
of the knot. In order to correct this problem, we used Kirkwood-Riseman 
formula Q to compute the friction coefficient 7 of every knot: since then 
v = qS/'y, we can find the specific time rescaling necessary to go from Monte 
Carlo time (tMc) to real time (tji): Ir = ( r yvMc/q£)tMC, where vmc is the 
velocity measured using tMC- Once we find the time conversion (one for 
every knot class) in the absence of the gel, then we apply it throughout our 
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simulations in the gel. 

Let us begin with the study of the dependence of the velocity with respect 
to the adimensional constant C = qSa/ksT for knot type 3i (if we assume 
that length fluctuations are very small, energy variations for the Metropolis 
algorithm can be expressed as integer multiples of fc^TC). We observe two 
distinct behaviors for the migration of knots (see fig. EJ) as a function of C. 

At high temperature or weak electric field, the distance of migration is 
linear as a function of C On the other hand, above a critical value of C cr n : 
the average speed of the knots is decreasing with C, in qualitative agreement 
with the experimental results. For our parameters, C cr u is located around 
0.4 for the 3i knot. Clearly, this value depends on the length of the knot 
and on its type (it can also depend on the gel parameter b if b < £, the 
typical size of the knot). If the electric field is strong, or temperature is 
low, knots tend to hang over obstacles and take a U-shape configuration. 
When a knot hits a gel rod, it can easily remain trapped there, because the 
probability of a backward step is very small, and it is a growing function of 
the temperature T. Trapping in U-shape conformations introduces plateaus 
in the migration distance as a function of time for individual knots, hence 
reducing the average migration velocity. On the same graph (fig. |3J), the 
drift velocity for a 81 knot is depicted as a function of C. The general 
behavior is similar to the 3i knot, but the drift speed at low C is higher than 
for the 3i knot (as it is the case in the experiments). But the most striking 
feature is that at C ~ 0.2 — 0.3 the two curves cross each other and the 3i 
is faster than the 81 for C > 0.3. This is also the observed behavior in the 
experiments. 

We investigate now the effect of the knot type for both weak and strong 
electric fields (C = 0.1 and C = 0.4 respectively) for eleven types of knots 
(3i, 4i, 5i, 52, 61, 62, 63, 7i, 72, 73 and 81) all consisting of 150 monomers. 
For each knot type, we extract the average velocity from the distance of 
migration vs. time curve. The velocity of migration is then plotted as 
a function of the measured ACN, that is related to the knot type. This 
plot is done first only for the high electric field case (C = 0.4) in fig. 0J 
We observe that there is a fairly linear relationship (except for the 3i knot) 
between the average velocity of knots and the ACN (measure of complexity). 
More complex knots migrate slower than simpler ones at strong electric fields 
(although much noisier than for weak fields) . These results are in agreement 
with experiments. A similar plot was done for the low electric field (C = 0.1) 
and the results are depicted in fig. |S] (filled squares). 

The intuitive view of a knot making its way through the pores of the 
gel would have as a consequence that more complex (thus more compact) 
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knots migrate faster than simple knots, because they are less disturbed by 
the gel. This is indeed what happens in weak electric fields. In the strong 
field regime our simulations are in agreement with experiments and show 
the opposite behavior, indicating that the knot collides with the gel and 
that the condition of self-avoidance makes the migration of compact knots 
around the gel strands more difficult. Somehow, parts of the knot have to 
go around other parts of itself, a process which is much longer for complex 
knots than for simple knots. 

The trapping of linear open polymers in U-shape conformations is ac- 
tually an artifact: indeed the slightest difference in the length of the two 
arms of the U gives rise to a net electric force that allows the polymer to 
slide around the obstacle. Simple local Monte Carlo moves do not allow 
to capture this dynamics. Yet, adding suitable non-local moves is enough 
to eliminate the artificious slowing down of the dynamics and cures the ex- 
ponentially long relaxation times of thermal activation around the obstacle 
[fil I28j. However, the closed knot topology of DNA in our numerical sim- 
ulations does not allow to introduce the long range moves. Moreover, we 
argue that exponentially long relaxation times present in our simulations 
affect, in first approximation, all the knots irrespective of their topology. 
In the ideal case of purely mechanical and frictionless unbinding, one can 
easily check that the knot complexity introduces at most a small logarith- 
mic correction according to which more complex knots would anyway unpin 
faster in contradiction with both experiments and our simulation. There- 
fore, the simulated absolute drift velocity v a f, s is affected by a time scale 
that is artificially stretched in essentially the same way independent on the 
knot class. So, by looking at the ratios of the absolute velocities, this time 
scale should in first approximation cancel out. In Figure [5] we plot the sim- 
ulated and experimental ratios as function of the ACN. In this 
Figure the open points represent data from Sogo et al.[2H) an d the filled 
symbols our simulation results. The ACN values for each knot were taken 
from Vologdskii et al. [201 ■ The agreement between experimental data and 
simulations is remarkable. In weak electric field, the transport properties of 
DNA knots are dominated by the hydrodynamics of the knots and the gel 
plays a minor role. At high electric fields, the knot-gel interaction is pre- 
dominant and it is responsible for the inversion of behavior. In particular, 
strong DNA-gel interactions enhance the effect of self-avoidance within the 
knotted DNA and self-avoidance must be included in simulations in order 
to reproduce the correct behavior. This is at variance with simulations of 
gel electrophoresis of linear DNA where self-avoidance is usually neglected 
because of the two following reasons: first, according to the repton model, 
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DNA crawls along tubes in the gel and it is in an elongated configuration 
where self-intersections are negligible. Secondly, it is often assumed that 
the drifting DNA can be considered to be in a semi-dilute solution regime, 
where polymers obey random walk statistics. Instead, in our case the con- 
servation of the knot topology during the dynamics imposes the inclusion of 
self-avoidance. 

We presented some results of a Monte Carlo simulation of DNA knots in 
a gel in three dimensions. In summary, our model allows to explain the high 
electric field behavior observed in experimental DNA knots gel electrophore- 
sis. The variation of the gel parameter does not change qualitatively the 
results. However, in a denser gel, the knots get stuck at shorter distances. 
In a more realistic modelization of the gel, the gel fibers should be allowed 
to break and let the knots migrate further. Varying the length of the knots 
would probably also bring some further insights into the mechanisms of the 
migration in a gel. 
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Figure Captions 



• Figure 1: Analysis by two-dimensional agarose gel electrophoresis of 
pH5.8 DNA (8749 bp) with the diagrammatic interpretation of the 
autoradiogram and the standard agarose electrophoresis (left panel). 
Reproduced from |2J with permission. 

• Figure 2: Elementary moves of the BFACF algorithm, a) creation of 
a handle (the opposite move destroys a handle) and b) corner flip. 

• Figure 3: Drift velocity in arbitrary units for the 3i (open circles) and 
8i (filled circles) knots as function of C = qSa/ksT. The lines are only 
guides for the eyes. 

• Figure 4: Linear relation between the electrophoretic drift velocity 
of the centre-of-mass as a function of the average crossing number 
(determined during the simulation) for different knots for low and high 
electric field (filled circles for C = 0.1, open circles for C = 0.4). 

• Figure 5: Distance of migration in gel for the experimental data or 
drift velocity for the simulated data vs. ACN (from [2U]) for knots 
from 3i to 8i. Open symbols are data from Sogo et al. |2j|, filled 
symbols are the simulated values of the drift velocity. Squares are 
for low electric field {£ = 0.6 V/cm and C = 0.1, respectively), while 
triangles are for high electric field (f = 5 V/cm and C = 0.4). The 
values were all normalized to their respective value for the knot 3i. 
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Figure 1: Analysis by two-dimensional agarose gel electrophoresis of pH5.8 
DNA (8749 bp) with the diagrammatic interpretation of the autoradiogram 
and the standard agarose electrophoresis (left panel). Reproduced from [2^ 
with permission. 



12 




Figure 2: Elementary moves of the BFACF algorithm, a) creation of a 
handle (the opposite move destroys a handle) and b) corner flip. 
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Figure 3: Drift velocity in arbitrary units for the 3i (open circles) and 81 
(filled circles) knots as function of C = qSa/ksT. The lines are only guides 
for the eyes. 
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Figure 4: Linear relation between the electrophoretic drift velocity of the 
centre-of-mass as a function of the average crossing number (determined 
during the simulation) for different knots for low and high electric field 
(filled circles for C = 0.1, open circles for C = 0.4). 
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Figure 5: Distance of migration in gel for the experimental data or drift 
velocity for the simulated data vs. ACN (from 20 ) for knots from 3i to 
81. Open symbols are data from Sogo et al. |2J, filled symbols are the 
simulated values of the drift velocity. Squares are for low electric field (E = 
0.6 V/cm and C = 0.1, respectively), while triangles are for high electric 
field (8 = 5 V/cm and C = 0.4). The values were all normalized to their 
respective value for the knot 3i. 
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